0^ 



X 



Multiphoton detachment of electrons from negative ions 

G. F. Gribakini and M. Yu. Kuchievi 
School of Physics, University of New South Wales, Sydney 2052, Australia 

Abstract 



(D . 

1^ ' A simple analytical solution for the problem of multiphoton detachment from 

m 

negative ions by a linearly polarized laser field is found. It is valid in the wide 
Q^' range of intensities and frequencies of the field, from the perturbation theory 

to the tunneling regime, and is applicable to the excess-photon as well as near- 
C^ ' threshold detachment. Practically, the formulae are valid when the number of 

O ■ photons is greater than two. They produce the total detachment rates, relative 

^~»' intensities of the excess-photon peaks, and photoelectron angular distributions 

for the hydrogen and halogen negative ions, in agreement with those obtained 
^.^ in other, more numerically involved calculations in both perturbative and non- 

(•~^ ■ perturbative regimes. Our approach explains the extreme sensitivity of the 

o' 

CN ■ multiphoton detachment probability to the asymptotic behaviour of the bound- 

O 



state wave function. Rapid oscillations in the angular dependence of the n- 
photon detachment probability are shown to arise due to interference of the 



C/3 

j>^' two classical trajectories which lead to the same final state after the electron 

O-c emerges at the opposite sides of the atom when the field is close to maximal. 
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I. INTRODUCTION 

In this paper we present an analytical solution to the problem of multiphoton detachment 
from a negative ion by a linearly polarized laser field. It gives very reliable quantitative results 
for a wide range of intensities and frequencies of the laser field, from the weak-field regime, 
where the process is described by the perturbation theory, to the strong fields where it proceeds 
as tunneling. The theory is valid when the number of photons n is large, but usually gives 
good results as soon as n > 2. We use it to calculate and examine various characteristics of the 
problem: the total multiphoton detachment rate, the n-photon detachment cross sections, the 
spectrum of excess-photon detachment (EPD) photoelectrons (the analogue of above-threshold 
ionization in atoms), and the peculiar photoelectron angular distributions. 

There are two important physical properties of the multiphoton detachment process: 
(i) The frequency of the laser field is much lower than the electron binding energy, 

^ < l^ol , (1) 

where Eq = —h?/2 is the energy of the bound state (atomic units are used throughout). 
This means that multiphoton detachment is an adiabatic problem. The external field varies 
slowly in comparison with the period of electron motion in the system. Therefore, the general 
adiabatic theory |jl|-^ is applicable. As long as the laser field is weaker than the atomic field, 
the detachment probability is exponentially small with respect to the adiabaticity parameter 
\Eq\/u ~ n. 

(ii) The process of multiphoton detachment takes place when the electron is far away from 
the atomic particle (see section |I|), at large distances 



r-R = J^^^^^, (2) 



Uy/lTj^ 



where 7 = ujk/F is the Keldysh parameter and F is the field strength. In the weak field 
regime, 7^1, Eq. (^ gives R ^ l/VcJ ~ K^^\/2n ^ 1, where k ~ 0.3 for a typical 
negative ion binding |ii^o| ~ 1 eV. In the strong field regime, 7 < 1, estimate (^ yields 



R ~ Jl /^ = k.~^JFo/F ^ 1, where Fq = k^ is the typical atomic electric field, F <^ Fq. 

The two features (i) and (ii) greatly simplify the multiphoton detachment problem. Owing 
to (ii), the final state of the electron can be described by the Volkov wave function ^ which 



takes into account the external field and neglects the atomic field. Moreover, the Volkov wave 
function describes explicitly the variation of the electron energy in the laser field. This makes 
it very convenient for application of the general adiabatic theory, as suggested by (i). 

Calculations based on the Volkov final-state wave function were first done by Keldysh f^. 
Subsequently, the idea was developed by Perelomov et el |^ and later reconsidered by Faisal 
and Reiss 0. This approach is usually supposed to give a correct quahtative picture of 
multiphoton processes. In this paper we demonstrate that, in fact, it produces very accurate 
quantitative results for the multiphoton detachment from negative ions. We re-examine and 
extend the Keldysh theory, paying particular attention to the following points. First, we show 
that the EPD can be accurately described by the theory. Originally, the theory was developed 
for low-energy photoelectrons |^ with kinetic energies much smaller than the binding energy. 
The present approach is valid at any photoelectron energy. Secondly, the angular distribution 
of photoelectrons is examined in detail. We show that a nontrivial oscillatory pattern of the 
angular distribution is caused by the simple and interesting physics. The photoelectron's 
escape from the atomic particle is most probable when the field reaches its maximum. There 
are two such instants in every period of the laser field T = 2tt/uj, say, t = and t = T/2. 
As a result, there are two classical trajectories which lead to the same final state of the 
photoelectron. Interference of the corresponding amplitudes gives rise to an oscillatory angular 
dependence of the detachment rate. There is a similar effect in the single-photon detachment 
in the presence of a static electric field, where the interference takes place between the two 
trajectories of the electron emitted up or down field [§|,|[. 

Estimate (0) leads to a further simplification of the problem, since the initial bound-state 
wave function of the atomic system should also be considered at large distances, where it can 
be replaced by its simple asymptotic form. The complicated behaviour of the wave function 
inside the atom, and the corresponding many-electron dynamics have little infiuence on the 
multiphoton detachment. In contrast, use of the wave function with incorrect asymptotic 
behaviour, e.g. that corresponding to the Hartree-Fock binding energy, introduces an error, 
which is exponentially large with respect to ^/n. Such sensitivity has been noticed in the 
perturbation theory calculations of the two- and three-photon detachment from H~ [^ . 

There has been a large number of papers where multiphoton detachment from the hy- 



drogen and halogen negative ions is investigated. Perturbation theory calculations include 
those based on the Hartree-Fock approximation |]rT|, adiabatic hyperspherical approach PUJ] , 
model potential [l^, a configuration-interaction procedure [Q, and the Lippmann-Schwinger 



equation |]I4| . There are also numerous non-perturbative methods, such as the Floquet close- 
coupling method [jl5|, complex- scaling generalized pseudospectral method [J16l, non-Hermitian 



Floquet Hamiltonian method jT^, and the i?-matrix Floquet theory [|l^,|T^. All the above 



methods rely on much more involved numerical calculations than those needed in our analyti- 
cal approach. However, we believe that the present theory provides accurate answers for most 
of the multiphoton detachment problems. For illustration purposes, we reproduce a variety of 
results obtained earlier, including the ra-photon cross sections, total detachment probability, 
EPD spectrum and photoelectron angular distributions for a large range of frequencies and 
intensities of the field (Sec. PTl] ). We believe that in some cases our results are more accurate 
than those obtained previously, due to the correct asymptotic behaviour of the bound-state 
wave function we use. 

The good accuracy we have achieved within the Keldysh-type theory is quite useful for the 
multiphoton detachment problem. On the other hand, its validity is very important for the 
development of an adiabatic theory of more complicated phenomena, such as double ionization 

The formulae obtained in this paper can be used to estimate probabilities of multiphoton 
ionization of neutral atoms. However, the infiuence of the Coulomb field of the positive ion 



on the wave function of the photoelectron cannot be neglected p3| , p^ , and our results for the 
multiphoton ionization would be less reliable. 

II. THEORY 

A. Basic equations 

Consider the removal of a valence electron from an atom or a negative ion by the laser field 
F(t) = F cos tut. The differential detachment rate can be written as the sum over n-photon 
processes [see Appendix]^, Eq. ( |A8|) ], 

dwn = 27r ^ \Apn\^S{Ep - Eq- nuj)——- , (3) 



where Apn is the amphtude of the n-photon process, 

Apn = ^ £ %{r,t)Vpmo{r,t)drdt , (4) 

\l/o(r, t) = e~*^°*$o(r) is the wave function of the initial electron state in the atomic potential 

f/(r), 

-,2 



y + Uir] 



<l>o(r) = ^o$o(r) , (5) 

Vpit) is the interaction with the laser field, 

V>(t) = -er • F(t) , (6) 

in the length gauge, e = — 1 for the electron, and \E'p(r, t) is the continuous spectrum solution 
of the time-dependent Shrodinger equation with the quasienergy Ep = p^/2 + F'^e^ /Aoj"^. It 
describes the outgoing photoelectron in the laser field with the translational momentum p, and 
F'^e^ j^LO^ is the electron quiver energy due to the field. The subscript n in A^^ reminds one 
that the amplitude must be calculated at Ep = Eq^ nuo provided by the energy conservation 
in Eq. (D. 

As we show below, the detachment probability is determined by the asymptotic behaviour 
of the bound-state wave function at large distances. This means that the role of electron corre- 
lations in the multiphoton detachment of a single electron is small, provided $o(r) represents 
correctly the asymptotic behaviour of the true many-electron wave function of the system, 

^7v(ri,...,rjv-i,r) ~ *^_i(ri, . . . , rAr_i)$o(r) , (7) 

where ^ n is the ground-state wave function of the A^-electron system, and \l'Ar_i is the wave 
function of the A^ — 1-electron atomic residue. 

If we neglect the infiuence of the atomic potential U{y) on the photoelectron, the final-state 
is given by the Volkov wave function. 



^p(r,t) = exp 



z(p + ki)-r-^y'(p + kiOW 



where k^ = e J*F(t')dt' is the classical electron momentum due to the field. By omitting the 
lower integration limit we mean that we set its contribution to zero, as if the integration is 
performed from — oo and the integrand is switched on adiabatically. For the Volkov function 



(I) this gives the same phase as in ||2T|], — |/o(p + iit')'^dt' + ipF/cu'^, and provides \I/p(r,t) 
with a convenient symmetry property with respect to inversion: 



^_p(r, t) = ^p(-r, t + T/2) exp(zEpT/2) . 



The wave function (^) satisfies the Shrodinger equation 



(9) 



' at 



y + V>(t) 



\1> 



(10) 



The neglect of the short-range potential U{r) for the photoelectron is justified in multiphoton 
processes, e.g., in the multiphoton detachment from negative ions, see end of Sec. |1IB| . 

Using the complex conjugates of Eqs. (§) and (p!OD, and id'^o/dt = Eq'^q, we transform 
amplitude (HI) into 



A 



1 r^ 



pn 



TJo 



E^ 



[p + ki 



<lo(p + ki)exp 



(p + ki/)^cit'-iEot 



dt , 



;ii) 



where $o(ci) is the Fourier transform of $o(r) 



*o(q) 



rfre-*'i'"$nfr) . 



(12) 



Note that in the velocity gauge, 

V>(t) = --A{t) . p+ ^A2(t) , A(t) = -cfY{t')dt' 
the Volkov wave function looks simpler. 



(13) 



*p(r, t) = exp 



ip-r-^J (p + ktOW 



(14) 



This gauge, which apparently 'leads to an analytical simplicity' 

.2\ ^T 



V = ^(^o-yl$o(p) 







exp 



(p + kt>ydt' - iEot 



dt 



(15) 



is though less physical than the length gauge in this problem. The amplitude (^ is not 
gauge invariant when U{r) is neglected for the final state (compare ([ll]) with (piSf )), except 
for the zero-range s-wave initial state, $o(i") = Ae''^^' /VAtt. The length gauge interaction 
(^) emphasizes large distances, where the bound-state wave function $o(r) has a well-defined 
asymptotic behaviour. We will see in the next section that this gives it a major advantage over 
the velocity gauge. In the limit u; ^ the length-gauge calculation reproduces the static-field 

result mm. 



B. Adiabatic approximation 

For multiphoton processes the integral over time in the amphtude ([Tl|), contains a rapidly 
oscillating exponent exp[iS{ujt)], where S{ujt) ~ 27m is the coordinate-independent part of 
the classical action 

S{ujt) = ^ J\p + kt^fdt' - Eot . (16) 

This makes the amplitude A-pn exponentially small, and the integral /g . . .dt should be cal- 
culated using the saddle-point method. The positions of the saddle points are given by 
dS{wt)/dt = 0, which yields 

{p + ktf + K^ = 0. (17) 

The saddle-point method in this problem has a simple and important physical contents. The 
two terms in the right-hand side of Eq. (p^ ) describe the energy of the electron in the initial 
and final states, Eq and (p + k()^/2, respectively. According to the general adiabatic theory 
0, the transition from the initial into the final state happens at the moment of time when 
their energies are equal. This is exactly the meaning of Eq. ([T7|). 



Note that condition (|T^ coincides with the positions of singularities of the Fourier trans- 



form $o(p + ki) in the amplitude (pHl) . Indeed, the general asymptotic form of $o(r) is 

$o(r) ~ Ar"-' exp{-tir)YiUr) , (18) 

r>l 

where u = Z/k, Z is the charge of the atomic residue {u = Z = ioi the negative ion), 
and r = r/r is the unit vector. It is easy to see that due to (|18|) the Fourier transform (p^) 
is singular at q^ = —k^. Using [27| we derive the following asymptotic form of $o(ci) for 



q — > ±iK 

(g^ + K 



$o(q) - A^A{±nuv) ^?5,^^1!+!^ ' (19) 



where (±)' = (±1)' corresponds to g -^ iz/t. 

Therefore, when the length-form amplitude is calculated by the saddle-point method, we 
do not need to know the behaviour of the bound-state wave function in the whole space. In 
contrast, when using the velocity-form amplitude (pISD , the value of the Fourier transform for 



the true final-state momentum p is needed. To calculate it one must know the exact wave 
function at all distances including r ~ 1. What makes the problem even harder, many-electron 
correlations become essential there. 

Equation (|l^) for the saddle-points presented explicitly as 

lp+—smutA +K^ = (20) 

defines complex values t^ where the transition from the bound state into the Volkov state takes 
place. Equation ( pO]) has two pairs of complex conjugate roots in the interval < Re {ut) < 27r. 
According to the general theory of adiabatic transitions p|, in the case when the final-state 
energy Ep is greater then the initial energy Eq, we should take into account the two saddle- 
points in the upper half-plane of complex t. 

Changing the integration variable to ut and substituting the asymptotic expression for the 
Fourier transform near the singularity (0), we can write amplitude (p!T|) as the sum over the 
two saddle points, 






iS{ujt)] d{uot) , (21) 



where the integral is taken over the vicinity of the /ith saddle point, p^ is the unit vector in 
the direction of p + (eF/cu) sintut^, and the two signs in (±) correspond to /i = 1, 2. Note that 
for the initial electron state bound by short-range forces, as in a negative ion, the integrand 
in (^) has no singularity (i/ = 0), and the application of standard saddle-point formulae is 



straightforward. Having the general case in mind, we will calculate the amplitude for arbitrary 
I/, taking into account the singularity at the saddle point. This is also useful if one wants to 
calculate the amplitude in the original form (^, without using the transformation which leads 
to Eq. (0). 

Using dS{ujt)/d{ujt) = [(p + k^)^ + k^]/2uj, we can re-write Eq. (^ as 

where = ut. In the vicinity of the saddle point 0^, 5"(0^) = 0, we have S"(0) ~ S'"(0^)(0 — 
(p^) . The contribution of this saddle point is then given by the following integral 

exp[z5(0)]^^^^ = -1— f ^^pMM1# , (23) 



[S'i(i>)Y ^"^ [S"MrJ (0-0^)^ 



which is calculated in Appendix 0. 

The explicit form of the action (|T^) is 

S{(j)) = n(j)- ^cos<i) - I sin 20 , (24) 

where z = e'^F^/Au^ is the mean quiver energy of the electron in the laser field in units of u, 
C, = eFp/uj"^ depends on the angle 6 between the photoelectron momentum p and the field F, 
and we put Ep — Eq = nuj due to the energy conservation in (^). Thus, we obtain the final 
expression for the amplitude by the saddle-point method, 






^"^ ' fj.=l,2 

where 



5 



(25) 



sin ut^ = (-e ± tJ8z{n - z) - e)/^z = s^, (26) 



coscjt^ = ±y^l - s2 = c^ , (27) 

S; = c,i^ + Azs,), (28) 

and the signs ± correspond to the two saddle points /i = 1,2. The usual definition of the 



spherical harmonics 28 










YU^,^) = - 


1 


_Qimifi/ 


-1)' 


n-\-\m 
2 



27r 



2[ + l (/-H)! l^^%H. ^^ (29) 

2 (/ + |m|)! ' ^ ' ^' ^ ^^ 



is generalized naturally to calculate Yim{Pf_i) for complex vectors by setting 



(p + kt) ■ F / pI 

COSt^ = ^ = \ 1 + ^r (30) 

(P + kt^F V K^ 



where the last equality is valid at the saddle points, p± is the component of p perpendicular 
to F, p±_ = psm6. The real physical angle 6 should not be confused with the complex angle 



■& from equations ( 129|) and (pOl). The azimuthal angle cp is the same in both cases. 

Using (P) and the symmetry of the spherical harmonics Y/m one can show that the ampli- 
tude (^) has the following symmetry properties: Apn -^ {—l)"'~^''Apn upon inversion p -^ — p 
{9^n — 9, ip—^ip + it), and Apn — > (— l)""^'"*"™ylp„, upon reflection in the plane perpen- 
dicular to the direction of the field {6 ^ n — 9). Consequently, the amplitude is zero for p 
perpendicular to the field, ii n + I + m is odd. 



It is easier to look at the physics behind Eqs. (P3|)-(PBD in the case when the photoelectron 
momentum is small, p <^ k. The following simpler formulae for the saddle points can be 
obtained from (|26|)-(p8D by setting ^ = 0, 



:F 



2 



smut^ = ±i7, cos utf, = ±\Jl + 7^, S'^ = i-f^l + 72—^ , (31) 

where 7 = kuj/F is the Keldysh parameter. Thus, for small photoelectron momenta the saddle 
points are tuti = i sinh~ 7 and ujt2 = n + i sinh~ 7, and the detachment takes place at the 
two instances when the external field is maximal, t = and T/2 on the real axis. Accordingly, 
the total amplitude (^Sf ) is the sum of the two contributions from these points. This results 
in oscillations in the photoelectron angular distribution, which we discuss in greater detail 
below. 

The original approach used in 0J^ was to expand Eqs. (p6D , (|27|) and the action (|^) 
in powers of p/n to the second order (see Sec. [IID2| ), thus obtaining corrections to (|3T1) . 
In this regime 7 remains the main parameter which determines the probability of multipho- 



ton detachment p9| . However, the applicability of the saddle-point result ( pSf ) is essentially 
narrowed by such expansion (Sec. p|) . 

The adiabatic nature of the problem allows us to estimate the radial distances which are 
important in the multiphoton detachment process. We have already seen that the saddle 
points in the integral in (0) coincide with the poles of the Fourier transform $o(q)- The form 
of $o(ci) at g ^ ±2K is given by the behaviour of $o(r) at r ^ 00. To estimate the essential 
distances look at Eq. (^21) . The range of (p where the integral is saturated is determined by 
|S"'(0^)(50)^| ~ 1, which gives 5(f) ~ |S"'(0^)|^^/^. The corresponding range of momenta p + kj 
is given by 

dq cos0^(50~ = — , (32) 

U \ 7 / it 



where we use Eqs. (pl|). The essential distances are obtained from r6q ^ 1, which yields 
estimate (|^). It is important that i? ^ 1 in both weak- and strong-field regimes. This makes 
the Keldysh approach valid for short-range potentials. There is another physical reason which 
helps to understand why the atomic potential can be neglected for the photoelectron. When a 
large number of photons is absorbed by the photoelectron, higher angular momentum partial 
waves are populated. The influence of the short-range potential upon them is small. For 



a given electron momentum p the important / values can be estimated as / ~ pR. In the 
perturbation theory regime this estimate yields / ~ {p/ K)^/n, which suggests that the spread 
of the probability to find the photoelectron with given / is described by a random walk of n 
steps. 

Estimate @) also explains the extreme sensitivity of the multiphoton detachment rates 
to the asymptotic behaviour of the bound-state wave function. Suppose a bound state wave 
function characterised by k' instead of the true k is used. The error in the amplitude @) 
introduced by replacing k by k' comes in as a factor exp[— A/ti?], where Ak = k' — k. The value 
of R is large, thus, even a small An can produce an exponential error in the amplitude. Using 
the perturbation-theory regime estimate of R we obtain the error factor of exp[— 2(AK/fi;) v^2n] 
for the detachment rate. 



C. Detachment rates 

The differential n-photon detachment rate for the electron in the initial state Im is obtained 
from Eqs. (|^) and ( pSj) after integration over (f and p, 



dwn 
dVL 



An \uj/ [l + \m\jl 



Pi 



\m\ 



1 + p'^ sin^ 9/k.'^ 






ey.Y>[-iCf,{i + zs^ 



(33) 



where p = J2{nuj — F'^/Aoo'^ + Eq) is the photoelectron momentum determined by the energy 
conservation, and z/ = for negative ions. According to the symmetry properties of Apn, the 
differential n-photon detachment rate is exactly zero at ^ = 7r/2 for odd n + I + m. 
The total n-photon detachment rate of the Im state is obtained by integrating (^), 



wi'"^^ = 2n r^smOdO 



/o a\l 

and if we are interested in the total detachment rate for a closed shell, the sum over m and 
the electron spin projections must be completed. 



I 

Wn = 2 ^ W. 
m=—l 



(Im) 



(35) 



The dominant contribution to this sum is given by the m = state, since it is extended along 



the direction of the field, see Sec. II D 2 



It is very easy to take the effect of fine-structure splitting into account. The two fine- 
structure components j = / ± ^ of a closed shell are characterized by different binding energies 
|ii^o| and values of k. The ra-photon detachment rate for the j sublevel is then given by 

dwj^'^ ^ 2j + l ^ dWn . . 

dn 21 + 1 ^, dn ' ^ ' 

m=—l 

which is exactly what one would expect from naive statistical considerations. 

Of course, one can easily obtain the total detachment rate by summing the n-photon rates 
over n. The smallest n is given by the integer part of [(|-Eo| + F^ / Au'^) / u] + 1. 

D. Limits 

There are two limits which can be usefully explored with the help of Eq. (^). The first 
is the perturbation theory limit, where the detachment rate is proportional to the nth power 
of the photon flux J = cF"^ / {Stiuj) , and the process is described by the generalized n-photon 
cross section 

dOn _ dWn n /oyN 

dQ~ dQ^ • ^^^^ 

The other is the low photoelectron energy limit studied earlier in [^0. It enables one to 
recover the static-field results p5| , p6| . 

1. Perturbation theory limit. 

To obtain the perturbation-theory limit, it is convenient to re- write the saddle-point equa- 
tion (^) for sincut^ in the following form: 



s, = j[p\\±i^k^+pI) , (38) 

where p\\ = pcosd is the momentum component parallel to the field. The weak-field regime 
7 3> 1 infers |s^| ^ 1, hence we obtain for cosu;t^ 



c, = ±Jl + sl^ -IS, + — + 0{s-') . (39) 



2s 



M 



Using (|38| ) and (^) to calculate the amplitude (^) for z/ = 0, and retaining only the leading 
term in c^ everywhere, except in c^ -|- zs^, where the second term is necessary, we arrive at 
the following n-photon detachment cross section 



dan pA^oj 

dn ~ 47r2v/2^ 
exp(pn/a;; 



(2/ 



(/-|m|)! 
'(/ + |m|)! 



P, 



m| 



1 + pi/K^ 



Tie 



ncu^ 



H? + p\ 



-l)"+'+"cosH 



(40) 



where p = y/2nij — k"^, c ~ 137 is the speed of hght, e = 2.71 . . ., and S is the momentum- 
dependent contribution to the relative phase of the two saddle-point terms in the amplitude, 

P\\ , p\\\/i^^+pI 



E = {2n + 1) tan"^ 



K^ +p5_ 



+ 



UJ 



(41) 



This phase varies with the ejection angle of the photoelectron from So = (2n+ 1) ta.n^^{p/K,) + 
pn/uj to —So, and can be quite large, even for the lowest n process, p ~ y/uj, Eq ~ y/n, thus 
producing oscillations in the photoelectron angular distribution. Note that in accordance with 
the general symmetry properties, the cross section is zero at 6 = tt/2, when n + I + m is odd. 



2. Low photoelectron energies and the static field limit. 

Another simplification of the general formula (^) is achieved when the energy of the 
photoelectron is low compared to the binding energy, p^ <^ h?. Then, following ||^,^] one can 
expand the action S{(f)^) and other quantities calculated at the saddle points in powers of p 
up to the 2nd order. For u = which corresponds to the multiphoton detachment from a 
negative ion, we obtain 
dwn pA'^uj-f 1 21 + 1 {I + \m\)\ 



dn 27r|Eo|Vl + 7^ (2l™l|m|!)2 4n {I - \m\)\ 



X exp < — 2- 



LU 



1 H sinh "'^7 

. 27 V 27 



exp 



sinh ^ 7 — 



7 



X exp 



7p2 sin^ 6\ ( p sin 6 



2\m\ r 



v/TT^ 



UJ 



l + (-l)"+'+'"cos 



\ p_ 

2Kpcos6y/l + 7^ 



wy 



(42) 



This formula coincides Eq. (53) of Ref. ||^. The cos(. . .) in the last square brackets of Eq. (^) 
appears due to the interference between the contributions of the two saddle points in amplitude 
(p5|), and is the analogue of cosS in Eq. (|40|). It determines the oscillatory behaviour of the 
angular dependence of the n-photon detachment rate, which would otherwise simply peak 
along the direction of the field, 6 = 0, or 6 = it, for m = 0. 



Formula (|42D also shows clearly that the detachment rate for the states with m 7^ is 
much smaller than that of m = 0, due to the factor {p sin 9 / ny^'"^^ . It comes from the leading 
term in the expansion of the associated Legendre polynomial P; (x) in Eq. ([33| ) at x ~ 1. 



As shown by Perelomov et al. |^, in the hniit u; — » 0, Eq. P^ ) allows one to recover the 
well-known formula for the ionization rate in the static electric field F [^, 

A^ (2/+l)(/ + |m|)! ^2Fo^^'-H-i / 2Fo\ 



for negative ion case z/ = 0. It has been shown recently [^ that the account of the polarization 



potential —ae^/2r^ acting between the outer electron and the atomic residue in the negative 
ion, changes the numerical pre-exponential factor in Eq. (^31). However, this correction is not 
very large, e.g., it increases the detachment rate for Ca^ by a factor of 2, in spite of the large 
polarizability a;(Ca) = 170 a.u. 

It is worth noting that the perturbation theory formula ( ^OD and the low electron energy 
limit (^2]) have a common range of applicability. If we use p -C k in the first, and take the 
perturbation theory limit 7 ^ 1 in the second, the two formulae yield identical results. 

III. NUMERICAL RESULTS AND DISCUSSION 

In this section we use the formulae we obtained within the adiabatic theory to calculate 
the photodetachment rates, EPD spectra and photoelectron angular distributions for H~ and 
halogen negative ions. These are so far the most studied species, which enables us to make 
comparisons with results of other calculations. Our aim is to show that our theory achieves 
good accuracy in describing multiphoton detachment in both perturbative and strong-field 
regimes. 

To apply the theory, all we need is the asymptotic parameters A and k of the corresponding 
bound-state wave functions. The values of A are tabulated in various sources, and we use 
those from WM. The values of k are calculated using the corresponding binding energies. 



K = J2\Eq\. They are taken from the electron affinity tables [^], or obtained by combining 
those with the fine-structure intervals of the atomic ground states ||3^, when we consider the 
detachment of pi/2 electrons from the halogens. 

In Fig. U we present the generalized n-photon detachment cross sections for H~ obtained 



by integrating the differential cross sections from Eq. (^OD with A = 0.75 and k = 0.2354 
over 6. The cross section has been multiplied by 2 to account for the two spin states [cf. 
Eq. (|36D with I = 0, j = 1/2]. The results of the perturbation theory calculations [|1^ are 



shown for comparison. In the latter the interaction of the electron with the atomic core was 
described by a model potential which accounted for the polarizational attraction between the 
electron and the atomic core, and was chosen to reproduce the binding energy of H~, as well 
as the electron-hydrogen scattering phaseshifts. Figure |l] shows that there is good agreement 
between our results and those of [J12l . We checked that even for n = 3 the difference does not 
exceed 20% at the cross section maximum. 

Laughlin and Chu note [jl2| that their model-potential results are close to those obtained in 
||T0| using the hyperspherical method which accounts for correlations between the two electrons 
in H^. They are also in agreement with the two-electron perturbation theory calculations of 



1^ and the recent i?- matrix Floquet theory calculations [0, which also take into account 



electron correlations. The main idea behind those approaches was to reproduce the negative 
ion wave function as correctly as possible at all distances, particularly near the atomic core. 
This idea is favoured by the experience gained in a number of problems, such as the single- 
photon detachment, electron-atom scattering, etc., where electron correlations are indeed very 
important. However, as shown above, the multiphoton problem under consideration proves to 
be different. Absorption of several quanta is dominated by large distances satisfying inequality 
(0). The complicated behaviour of the wave function inside the atomic core turns out to be 
inessential. This is the main reason for the good agreement we observe in Fig. |I|. 

To check our theory in the non-perturbative regime, where one must use Eq. (|33D, the 
EPD spectra of H~ for the three large field intensities / = 10^*^, 5 x 10^°, and 10^^ W/cm^, 
of the 10.6-/im radiation, uo = 0.0043 a.u., are presented in Table |. For these parameters the 
electron quiver energy, or the ponderomotive energy shift, in units oi u, z = F'^/Au^ = 0.894, 
4.472, and 8.945, and the Keldysh parameter 7 = 1.895, 0.847, and 0.599, respectively. For 
given uj absorption of a minimum of 7 photons is required. The ponderomotive threshold shift 
changes this number to nmin = 8, 11, and 16. The calculation of the detachment rates from 
Eqs. (|33D -(^6|) has been done using Mathematica ||3l]]. For the smallest intensity the lowest 
EPD peak n = 8 dominates the total detachment rate, whereas for the higher intensities many 
peaks in the EPD spectrum can be observed. 

The detachment rates in Table | are compared with those obtained in the non-perturbative 
calculations of Telnov and Chu [1^. They describe their method as a complex-scaling general- 



ized pseudospectral technique applied to the solution of the time-independent non-Hermitian 
Floquet Hamiltonian for the complex quasienergies, and use the accurate model potential 
from [J12l to describe the interaction of the electron with the atomic residue. 

There is a good overall agreement between the two calculations. The discrepancy usually 
does not exceed a few per cent, and is slightly larger for higher EPD peaks and smaller 
field intensities. The latter is somewhat puzzling, since there is a good agreement in the 
perturbation-theory limit for the 7-photon cross section at a; = 0.0043 a.u.: 

0-7 = 3.537 X 10"^°° cm^S^ (Eq. (40) integrated over angles), 
ay = 3.639 X 10-^°° cm^S^ (result of W%). 



In Fig. ^ we show the angular dependence of the photoelectron peaks for n = 16, 17, 18, 
and 19, at / = 10^^ W/cm^. We have checked that their shapes, as well as those for other n 
and intensities, are practically identical to the angular distributions presented in Figs. 5-7 of 
P!7| . Also shown in Fig. ^ are the differential detachment rates obtained from Eq. (^^. It 
works quite well for two lowest n, but the agreement becomes poor with the increase of the 
photoelectron energy, e.g., for n = 19, where p/n ^ 0.75. 

It is worth stressing again that the remarkable oscillatory behaviour is caused by the inter- 
ference of the two saddle-point contributions in Eq. (|33|) , or in other words, the interference 
between the electron waves emitted at the two instants separated by T/2, when the field 
reaches its maximum. The geometrical phase difference which determines the oscillations of 
cos(. . .) in Eq. ( ^2]) can be calculated classically. Suppose that the electron is considered free 
at the moment when it escapes the atomic particle. Its classical coordinate is then given by 
r(t) = j^'ktidt' = (F/cu^) cos cut. At the two instants t^ when the adiabatic transition takes 
place we have 



where Eq. (|3l[) is used for small momenta p <^ k. Note that though t^ are complex, the 
corresponding electron coordinates are real. These points located at the opposite sides of the 
atomic particle, are sources of the two electron waves emitted at the angle 6 with respect to 
F. The geometrical phase is obtained by multiplying the base |r(ti) — r(t2)| by the projection 
of the electron momentum on the the direction of the field pcos6. 



Our results for halogen negative ions are presented in Fig. ^ and Table ^ They have been 
obtained from Eqs. ppj), ( PU] ) for comparison with the perturbation theory calculations [|rT| 
at the Nd:YAG laser frequency uj = 0.0428 a.u. In that work the non-relativistic Hartree-Fock 
wave functions of the valence np electrons were used, together with experimental threshold 
energies. The photoelectron was described in the plane wave approximation. This approxi- 
mation is equivalent to our use of the Volkov wave function in the perturbation-theory limit. 



As shown in the earlier works by Crance [Q, the multiphoton detachment results obtained in 
the plane wave approximation are close to those obtained using the frozen core Hartree-Fock 
wave functions of the photoelectron. 

The shapes of angular distributions presented in Fig. |^ are quite close to those in Fig. 2 
of Ref. |Tl|, although quantitative comparison is not feasible due to the use of an arbitrary 



vertical scale in Wl 



The absolute values of the n-photon detachment cross sections in from our calculations and 
TT| compare reasonably on a logarithmic scale for all cases shown in Table |l[ However, there 



is a systematic discrepancy. To find its origin let us recall that the multiphoton detachment 
rate is very sensitive to the asymptotic behaviour of the bound-state wave function (see end 
of Sec. [II Bl) . In [^ the Hartree-Fock wave functions have been used. Their asymptotic 



behaviour exp(— /thf^) is different from the correct exp(— /tr), based on the experimental value 
of K. Thus, to account for the discrepancy in Table ||, the Hartree-Fock based results should 
be multiplied by the factor 

~ exp[2(fi:hf - f^)R] , (44) 

where, according to (^, i? ~ l/^/cJ- Formula ( PD shows that when /thf > k, the Hartree-Fock 
based calculations underestimate the detachment rate, while for /thf < ^ they overestimate it. 
The Hartree-Fock values of /thf are 0.602, 0.545, 0.528, and 0.508, for the outer np subshell 
of F~, CI", Br", and I~, respectively. Examination of the lowest n cross sections throughout 
Table |l| shows that the qualitative explanation of the discrepancy based on (|4^) is correct. 
For example, for F~ where k^ = 0.6 and n = 0.5, formula (^4]) gives 2.6, whereas the ratio 
of the 3-photon detachment cross sections for F~ , j = 3/2, in Table || is 4.3. Also, the best 
agreement in Table |I| is achieved for Br~, j = 1/2, where km is very close to the correct 
value. Therefore, we conclude that the incorrect asymptotic behaviour of the Hartree-Fock 



wave functions can produce significant errors in the multiphoton detachment rates. This 
must be kept in mind when comparisons are made between different n-photon detachment 
calculations |jl9l . 



IV. SUMMARY 

The main result of our work is that the adiabatic theory approach to the multiphoton 
problems originally suggested by Keldysh, is more powerful and accurate than is generally 
believed. It yields accurate multiphoton detachment rates for negative ions, and reveals a 
number of interesting details about the physics of the problem: the role of large distances and 
asymptotic behaviour of the bound-state wave function, and the origin of oscillations in the 
angular distribution of photoelectrons. The formulae obtained in the paper allow one to make 
simple and reliable estimates of the n-photon detachment rates for n > 2 in both perturbative 
and non-perturbative regimes. 
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APPENDIX A: CALCULATION OF TRANSITION RATES IN A STRONG 

PERIODIC FIELD 

Suppose the system is in the initial state 

of the time- independent Hamiltonian Hq, and a periodic field V{t) = V{t + T) is turned on 
adiabatically. We assume that this field can be strong, so that the lowest-order perturbation 
theory is inapplicable. The time-dependent wave function of the system 

^^ = [Ho + V{m (Al) 



can be presented as the sum 

^(t)=V^o(t) + E«A(t)^A(t) (A2) 

A 

over the set of eigenstates V'a(^) of the total Hamiltonian, 

which represent the possible final states of the system, a\{t) being the amplitude of finding 
the system in one of these states. In Eq. ( [A^ ) we assume that axit) — >• at t — ;> — cx), and 
the rate of the transition ipQ -^ ipx is given by d\ax{t)\'^/dt. 

According to the Floquet theorem each state ipxit) = e~*^^*0A(^) is characterized by its 
quasienergy Ex and the corresponding periodic quasienergy wave function (pxit) = 4>\{t + T), 
found from 

^^ = [Ho + V{t) - Ex]<Px . 

At any given t the quasienergy wave functions form a complete orthonormal set, {ipxl'ipx') = 
{4>x\4>x') = Sxy- 

After inserting \l'(t) ( |A^ ) into Eq. (|A1|) and projecting it onto the state {ipxit)\, we arrive 
at 

day 



A - _,7W,, ('y-^lT/mU/,„m\ - _^JExt^-iEot 



z{Mt)\V{t)\Mt)) = -^e*^^*e-^«*(0A(t)|K(t)|0o) . (A3) 



dt 
The last matrix element is a periodic function of time 



,(t)|F(t)|0o)=Ee"*""*^An, (A4) 



n 



where u = 2tt/T, and 



Axn = Tf [ (0A(t)|K(t)|0o)e^""*rft . (A5) 

I JO 



Using ( |A4| ) we re-write ( [A 3D as 



dt V 



and find 



ax = —r- = — > — A 

J dt '^ Ex — Eq — nw — IT) 



Xn 



where the energies E\ have been given an infinitesimal shift —ir] to make J^ . . .dt converge at 
— cxD. The probabihty is given by 

IoaI^ = V -7-= ^-^^^-- — r + oscillating terms , 

and the rate is 

d. ,,,2 V- 2r7e2''* 



^a(^)p = e .. :'\.....m - 



|2 



dV ^^ '^ ^ {E\ — Eq — nwY + rf^ 

where we dropped the oscillating terms since they do not contribute to the transition rate after 
we average it over a period. Finally, we take the limit r^ ^ using the following representation 
of the (5-function, 

lim -^ = 2ti5{x) 

and obtain 

j^\a^{t)\^ = 2TiY.\A^^\H{E^-E^-nw) , (A6) 

where the amplitude A\n given by Eq. ( |7^ ) can be written as 

-j rp 

Axn = j;l {Mt)\V{t)\Mt))dt , (A7) 

due to the energy conservation Ex — Eq = nw implied by the ^-function. This amplitude 
describes the n-quantum process, and the total transition rate ( |A^ ) is the sum over all such 
processes. If the spectrum of A is continuous, the differential transition rate dw\ is proportional 
to the corresponding density of states, 

dwx = 27iJ2 \Axn\^S{Ex -Eo- nw)dpx . (A8) 

n 

APPENDIX B: SADDLE-POINT METHOD FOR INTEGRALS WITH A 

SINGULARITY 

Consider the integral 

J= [ g{x)exp[-\f{x)]dx (Bl) 

c 



for A — i> oo. In this case it is well known ||3^ that the integration contour C should be 
deformed to go through the saddle point xq where f'{xo) = 0. The vicinity of this point gives 
the main contribution to the integral. If the function g{x) is not singular at x = Xq, the 
integral (pl| ) is evaluated as 



J^9{xo)^jjj^^exp[-Xf{xo)] . (B2) 

If g{x) has a singularity at x = xq, e.g. g{x) = (x — xq)^'^, the saddle-point answer has to 
be modified. Consider the following integral, 



By using the transformation |2^ 

1 



(x - xoY T{u) 
we turn ( p3|) into the double integral 



1 r°° 

— d^C'expl-^x-xo)], (B4) 

[ly) Jo 



j^ d^^-' J exp[-A/(x) - e(x - Xo)]dx . (B5) 



Calculating J . . .dx by means of (|B2|) and then integrating over ^ we obtain for A — i> oo, 



•'""'"iff Vv5"'^"<^°""''^="^'-'^<'^">' ■ '^'' 



For z/ — > we, of course, recover (|B2|) . 
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FIGURES 
FIG. 1. Frequency dependence of the generalized n-photon detachment cross sections for H~, 
n = 5, 6, 7. Sohd curve: present calculation, Eq. ( ^0| ) integrated over angles; open circles: pertur- 
bation theory calculations of Laughlin and Chu ll^ . 



FIG. 2. Differential n-photon detachment rates of H in the strong laser field, uj = 0.0043 a.u., 
/ = 10" W/cm2, z = 8.945, 7 = 0.599. Sohd curve: "exact" saddle-point calculation, Eq. (|3|); 
dashed curve: low photoelectron energy limit, Eq. (^). Channels with n < 16 are closed. 

FIG. 3. Differential n-photon cross sections for the electron detachment from the halogen negative 
ions, which leaves the atom in the '^Ps/2 or '^Pi/2 states, Eqs. (4C), ([3^). 



TABLES 



TABLE I. The EPD spectra of H^ in the strong laser field of a; = 0.0043 a.u. The detachment 
rates calculated by our saddle-point method (SP), Eq. (p3|), A = 0.75 and k = 0.235, are compared 
with the non-perturbative results by Telnov and Chu [0]. 
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TABLE II. Comparison of the n-photon detachment cross sections from the halogen negative ions 
obtained by the saddle-point method (SP), Eqs. (pO|), (p6|), with the perturbation theory calculations 
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by Crance [|ll| at u; = 0.0428 a.u. For each n, log (Tn is shown, ali being in units of cm s 
j = 3/2 and 1/2 for the ^i-*3/2 and ^Pi /2 final states of the atom. 
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Figure 2. 
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